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We propose a new Dark Energy parametrization based on the dynamics of a scalar field. We use 
an equation of state w = (x — l)/(x+ 1), with x = Ek/V, the ratio of kinetic energy _Efe = 4> 2 /2 and 
potential V. The eq. of motion gives x = (L/6)(V73.ff 2 ) and with a solution x = ([1 + 2L/3(1 + 
y)] 1/2 - 1)(1 + j/)/2 where y = p m /V and L = (V'/V) 2 {1 + q) 2 , q = 4>/V . Since the universe 
is accelerating at present time we use the slow roll approximation in which case we have \q\ <C 1 
and L ~ (V' /V) 2 . However, the derivation of L is exact and has no approximation. By choosing 
an appropriate ansatz for L we obtain a wide class of behavior for the evolution of Dark Energy 
without the need to specify the potential V. In fact w can either grow and later decrease, or other 
way around, as a function of redshift and it is constraint between — 1 < w < 1 as for any canonical 
scalar field with only gravitational interaction. Furthermore, we also calculate the perturbations of 
DE and since the evolution of DE is motivated by the dynamics of a scalar field the homogenous 
and its perturbations can be used to determine the form of the potential and the nature of Dark 
Energy. Since our parametrization is on L we can easily connect it with the scalar potential V(4>). 



I. INTRODUCTION 



In the last years the study of our universe has received 
a great deal of attention since on the one hand fundamen- 
tal theoretical questions remain unanswered and on the 
other hand we have now the opportunity to measure the 
cosmological parameters with an extraordinary precision. 
Existing observational experiments involve measurement 
on CMB [1] or large scale structure LSS [2] or supernovae 
SNla [3], and new proposals arc carried out [4]. It has 
been established that our universe is fiat and dominated 
at present time by Dark Energy "DE" and Dark Mat- 
ter "DM" with Q DE ^ 0.73, Cl DM ^ 0.27 and curvature 
fij. ~ —0.012 [1]. However, the nature and dynamics 
of Dark Energy is a topic of mayor interest in the field 
[5]. The equation of state "EOS" of DE is at present 
time w ~ —0.93 but we do not have a precise measure- 
ments of w(z) as a function of redshift z [1]. Since the 
properties of Dark Energy are still under investigation, 
different DE parametrization have been proposed to help 
discern on the dynamics of DE [6] -[12]. Some of these DE 
parametrization have the advantage of having a reduced 
number of parameters, but they may lack of a physical 
motivation and may also be too restrictive. Furthermore, 
the evolution of DE may not be enough to distinguish be- 
tween different models and the perturbations of DE may 
be fundamental to differentiate them. 

Perhaps the best physically motivated candidates for 
Dark Energy are scalar fields which can be minimally 
coupled, only via gravity, to other fluids [11, 12] or can in- 
teract weakly in interacting Dark Energy "IDE" [15, 16]. 
Scalar fields have been widely studied in the literature 
[11, 12] and special interest was devoted to tracker fields 
[11] since in this case the behavior of the scalar field <f> is 
very weakly dependent on the initial conditions at a very 
early epoch and well before matter-radiation equality. In 
this class of models the fundamental question of why DE 
is relevant now, also called the coincidence problem, can 
be ameliorated by the insensitivity of the late time dy- 



namics on the initial conditions of 4>. However, at present 
time we are not concern at this stage with the initial con- 
ditions but we work from present redshift z = to larger 
values of z in the region where DE and its perturbations 
are relevant. In this case the conditions for tracker fields 
do not necessarily apply. Interesting models for DE and 
DM have been proposed using gauge groups, similar to 
QCD in particle physics, and have been studied to un- 
derstand the nature of Dark Energy [13] and also Dark 
Matter [14]. 



Here we propose a new DE parametrization based on 
scalar fields dynamics, but the parametrization of w can 
be used without the connection to scalar fields. This 
parametrization has a reach structure that allows w to 
have different evolutions and it may grow and later de- 
crease or other way around. We also determine the per- 
turbations of DE which together with the evolution of 
the homogenous part can single out the nature of DE. 
With the underlying connection between the evolution 
of w and the dynamics of scalar field we could determine 
the potential V ((/))■ The same motivation of parame- 
terizing the evolution of scalar field was presented in an 
interesting paper [10]. We share the same motivation but 
wc follow a different path which has the same number of 
parameters but it has a richer structure and it is easier 
to extract information on the scalar potential V(4>). 



We organized the work as follows: in Sec. I A we give 
a brief overview of our DE parametrization. In Sec. II 
we present the dynamics of a scalar field and the set 
up for our DE parametrization presented in Sec. III. We 
calculate the DE perturbations in Sec. IV and finally we 
conclude in Scc.V. 
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A. Overview 

We present here an overview of our w parametrization. 
The EOS is 



V 

w = — 
P 



x - 1 
x + l 



(1) 



with x = Ek/V the ratio of kinetic energy and poten- 
tial. The equation of motion of the scalar field gives (c.f. 
cq.(17)), 



1 



2L 



(2) 



where L = (V'/V) 2 A, y = p m /V the ratio of matter 
and V and A = (1 + q) 2 ,q = 4>/V. Eq.(2) is an exact 
equation and is valid for any fluid evolution and/or for 
arbitrary potentials V{4>). 

The aim of our proposed parametrization for L, y is 
to cover a wide range od DE behavior. Of course other 
interesting paramctcrizations are possible. From the dy- 
namics of scalar fields we know that the evolution of w 
close to present time is very model dependent. For exam- 
ple, in the case of V = V (f)~ 2 / 3 , used as a model of DE 
derived from gauge theory [13], the shape of w(z) close to 
present time depends on the initial conditions and it may 
grow or decrease as a function of redshift z. Of course if 
we change the potential the sensitivity on the choice of V 
and initial conditions will vary a lot. We also know that 
tracker fields are attractor solutions but in most cases 
they do not give a negative enough w [11]. The dy- 
namics of scalar fields with a single potential term, for a 
wide class of models, gives an accelerating universe only 
if A = V'/V -)• or to a constant |A| with w = -l + A/3 
[12]. In this class of models the EOS, regardless of its ini- 
tial value, goes to a period of kinetic domination where 
w ~ 1 and later has a steep transition to w ~ -1, which 
may be close to present time, and finally it grows to w in 
a very model and initial condition dependent. Further- 
more, if instead of having a single potential term we have 
two competing terms close to present time, the evolution 
of w(z) would even be more complicated. Therefore, in- 
stead of deriving the potential V from theoretical models 
as in [13] we propose to use an ansatz for the functions 
L, y which on the one hand should cover as wide as pos- 
sible the different classes of DE behavior with the least 
number of parameters, without sacrificing generality, and 
on the other hand we like to have the ansatz as close as 
possible to the know scalar field dynamics. We believe 
that using or model will greatly simplify the extraction 
od DE from the future observational data. We propose 
therefore the ansatz (c.f. eq.(41)) 



L =L a + L iy tf(a) = L Q + L lV l 



1 + (a/a t ) k 



/(«) 



l + {a/a t ) k l + [(l + z t )/(l + z)] fc 



(3) 
(4) 



where L a , L\ are free parameters giving w and W\ = 
w(z 3> z t ) at early times, f(z) is a function that goes 
from/(z = 0) = l/(l+(l+z t ) fe ) at z = to f(z > 1) = 1 
and Zt sets the transition redshift between w and w\ (a 
subscript o represents present time quantities) while k 
the steepness of the transition and £ takes only two val- 
ues £ = 1 or £ = 0. We show that a steep transition of w 
has a bump in the adiabatic sound speed c 2 which could 
be detected in large scale structure. Since the universe 
is accelerating at present time we may take the slow roll 
approximation where \q\ -C 1,A ~ 1 and L ~ (V'/V) 2 . 
However, the derivation of L in eq.(2) is exact and has no 
approximation. We will show in section III that w can 
have a wide range of behavior and in particular it can de- 
crease and later increasing as a function of redshift and 
viceversa, i.c the shape and steepness are not predeter- 
mined by the choice of parametrization. Of course we 
could use other parametrization since the evolution of x 
and w in eqs.(l) and (2) are fully valid. There is also no 
need to have any reference to the underlying scalar field 
dynamics, i.e. it is not constraint to scalar field dynam- 
ics. However, it is when we interpret x = (f> 2 /2V with 
L = (V' jV) 2 A and the ratio y = p m /V that we connect 
the evolution of w to the scalar potential V ((/>). 



II. SCALAR FIELD DYNAMICS 

We are interested in obtaining a new DE parametriza- 
tion inferred from scalar fields. Since it is derived from 
the dynamics of a scalar field <j> we can also determine its 
perturbations which are relevant in large scale structure 
formation. We start from the equation of motion for a 
canonical scalar field cf>(t, x) with a potential V(<p) in a 
FRW metric. The homogenous part of </> has an equation 
of motion 



cf> + 3H<j> + V'(<(>) =0 



(5) 



where V' = dV/dcj), H = a/a is the Hubble constant, a 
is the scale factor with a /a = 1 + z and a dot represents 
derivative with respect to time t. Since we are interested 
in the epoch for small z we only need to consider matter 
and DE and we have 



ZH 2 = Pn 



P4> 



(6) 



in natural units 8irG = 1. The energy density p and 
pressure p for the scalar field are 



P0 = \j? + V, P ^ = \^ 2 - V 

and the equation of state parameter "EOS" is 

VjL = 2 /2 - V = x-l 
W ~ P4> 4> 2 /2 + V x + l 

where we have defined the ratio of kinetic energy and 
potential energy as 



(7) 



(8) 



2V 



(9) 
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The value of x gives w or inverting eq.(8) we have x = 
(w+l)/(w — 1). The evolution of x(z) > determines the 
Dark Energy w in the range — 1 < w < 1. For growing x 
the EOS w becomes larger and at x S> 1 one has w ~ 1 
while a decreasing x has w approaching -1 for x = 0. 
In terms of x and y = p m /V~we have 



V(x + 1), p^ = V(x-l), = Vl/ (10) 



and 



1 + x 
1 + x + y 



y 



1 + x + y 



3H 2 = Pm + Pc/l = V(l + a: + y) 
We can write 



(11) 
(12) 



2V 'SH 2 6V 2 



V V' 2 9 V L 

< 1 + «»=3lP6 < U » 



using 



V + V'(l + q) 



3H 



3H 



and we defined 



V 



:[-) (l + g) 2 = A 2 A 



(14) 



(15) 



x- V ' 
X=y, q 



V 



A = (i + q y 



with q > —1. Since the r.h.s. of cq.(13) still depends on 
x through H we use eq.(12) and eq.(13) becomes then 



L 



6(l + x + y) 
which has a simple solution 



(16) 



(17) 



Eq.(17) sets our DE parametrization as a function of L 
and y. For small q and A is close to one and the quan- 
tity L gives direct information on the potential and its 
derivative. 



1. Dynamical evolution of x and y 

Differentiating x and y = p m /V w.r.t. time we get the 
evolution 

* (t-xV) =^-{q-x) (18) 



V V J V 

x — q 



6Hx 



1 + 9 



(19) 



and 



Pm yV'4> (2x-q-l\ 

y = — - -ttt = iH v \ — i i ) ( 20 ) 



v v 2 



1 + 5 



were we used eqs.(13), (14) and p m = —3Hp m (l + w rn ) 
with w rn = 0. The dynamical system of a scalar field with 
arbitrary potential V was studied in [12], and the critical 
points with x = in eq.(19) and constant x are are = 
and cf>— xV' — or equivalcntly x = 0, x = q, respectively 
The first case, = 0, implies x = 0,w = — 1 and a 
constant V(4>) with f2^ — > 1. At the same time, eq.(20) 
gives y = —3Hy with has a solution y = yi(a/ai)^ 3 — > 
and f2„,, — > 0. In the second case, q = 0/V = x, 
depending on the value of q the quantities x — q and w 
will take different constant values, and for w < w rn = 
(i.e. x < 1) we will have an increasing in — > 1 [12]. 
Setting q = x in eq.(20) we get ?) = 3Hyw with w = 
(x — 1)/(.t + 1) constant giving a solution 



y = Vo\ — 

. a„ 



'Atv 



(21) 



The critical points y = in eq.(20) arc y — and g = 
2:r — 1 with y constant. In the first case we have Cl m = 
and fl^ — 1 while in the second case eq.(19) becomes 
x = 3H(l — x)/x with a solution 



\l-x) 



e x (l-x ) - 



(22) 



At large values of a the l.h.s. of eq.(22) vanishes, x — > 
l,q = 2x — 1— >l,ui— >0 and f2 TO , f2^ are also constant 
with = 1 — O m (for a generic barotropic fluid the 
critical point would have been w — > w m which in our 
case is w m = 0). Clearly if x = q ar x = we are at 
a critical point x — in cq.(19) with constant x,w but 
for x ^ q,0 the system evolves and w is in general not 
constant. 

However, in the present work we do not want to study 
the critical points but the evolution of x close to present 
time when the universe is accelerating with x close to 
zero (w close to -1) but not exactly zero with 0^0 and 
<j) ^ xV' . We can assume that V — V'4> < 0, since we 
expect <f> to roll down the minimum of its potential, and 
the second term in eq.(18) is then negative, while <f> can 
take either sign. In the region where 4>/V < x we have 
x > while for 4>/V > x we have x < 0. If we take what 
we call a full slow roll defined by <j> = and 3H(f> = — V' 
then eq.(18) becomes 



x = -(1 + z)Hx z = QHx 2 



(23) 



which is positive definite, i.e. x > 0,x z = dx/dz < 0. 
The solution to eq.(23) is 



x(a) 



1 + 6x Ln(a /a) 1 + 6x Ln(l + z) 



(24) 



Therefore if the condition <fi = Oor \q = <j)/V'\ <C x 
is satisfied eq.(24) gives a decreasing function for x as a 
function of z and therefore w(z) also decreases. However, 
we do not expect to be in a full slow roll regime and when 
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x is small, e.g. w < —0.9 one has x < 0.05, the slow roll 
condition |</>| < \V'\ does not imply that <j> <C xV and 
the sign of x can be positive or negative depending on 
the sign and size of q = 4>/V compared to x and x can 
either grow or decrease. The value of q parameterizes the 
amount of slow roll of the potential and a full slow roll 
has q — but we expect to be only in an approximate 
slow roll regime with \q\ <C 1 and A ~ 1. We will discuss 
further the value of q in section IV. 



A. Evolution 

Taking the differential of eq.(16) we have 
—x dy dL 



dx = x y dy + x^dL 



l + 2x + y 6(1 + 2x + y) 



(25) 



where the subscript represents a derivative, e.g. x y = 
dx/dy, or in terms of derivative w.r.t. the redshift z we 
have 



X z — x yVz + XlL z 



with 



6(1 + 2a; + y) 
1 



XL = 



l + 2x + y 



(26) 

(27) 
(28) 



Clearly xl is positive definite while x y < 0. The deriva- 
tive of w is 



UK 



W X X Z 



with 



2x 



> 0. 



(29) 



(30) 



In general we can assume that DE rcdshifts slower than 
matter, at least for small z, since DE has w < and 
matter w m = 0, so y = p m /V is a growing function of 
z, i.e. y z > 0. For L constant we see from eq.(26) that 
x z = x y y z < and x will decrease and so will w. On 
the other hand for y constant we have x z = x y dL and an 
increase on L gives a larger x and w. 

We present in the appendix the dynamical equations 
of L and q and the limits of x from eq.(17) the following: 
in the limit L <C 1 we have from eq.(16) with 



L 



-1 



L 



3(1+ tf) 



(31) 



and therefore x — > 0, w — > — 1 as L — > 0. For y ^> 1, with 
L constant, we have 



6y' 



w 



-1 



(32) 



and again we have x — > 0, w 
with L constant, we have 



- 1 as y — > oo . For y <C 1 , 



L 1 L 

6' W = - 1+ 3 



(33) 



giving a constant x and w. For L< 1, with y constant, 
we have 



L 

6' W 



1 



(34) 



giving a constant x and iu. Finally, the limit L/y — > L\ 
constant with y 3> 1 has a constant a; and w with 



6 : 



Li -6 
u> = = — 1 



Li+6 



2Li 
6 + Li' 



(35) 



As we see from eqs.(32), (33) and (35) all these limits are 
obtained from eq.(31), i.e. eq.(31) is then valid for the 
limits L < 1 or y > 1 or y < 1. Clearly depending on 
the choice of L we can have a decreasing or increasing 
i,«iasa function of redshift. 

As we see from eq.(32) for y ^> 1 (i.e p m ^> V) valid 
at larger z we can estimate the evolution of L/y from 
cq.(20) and (A7) giving 



d(L/y) 3HL I H 



1 



and using 
H 



3H 2 



dt 



3H 2 1 



2x 



2(1 + x + y) 



(36) 



(37) 



we have d{L/y)/dt > for q > (1 + x(l + 2y + 2x))/(l + 
2y + 3x). 



1. Late time attractor Solution 

The evolution of scalar field has been studied in [12] 
and the late time attractor for scalar fields leading to 
an accelerating universe with — > 1 requires w < 
-1/3, L < 2. In the limit |A| > 1 one has [12] with 
Pm < P<p and 



6H 2 



L 
6' 



V 
3#2 



L 
6" 



giving 



2V 



L 



V 



L 

- = -1+3 



(38) 



(39) 



If we take the limit 5 C 1 and in eq.(31) we recover w 
in eq(39). However, for large z we expect y to increase 
and eq.(39) would not longer be valid. In this region 
we should use eq.(31) or the full value x from eq.(17) in 
cq.(8). 
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III. SCALAR FIELD DE PARAMETRIZATION 

In order to have an explicit parametrization of Dark 
Energy we need to either choose a potential V(<j)) or take 
a parametrization for L and y. Of course if we want to 
study a specific potential V we would solve the equation 
of motion in eq.(5). However, the aim here is to test a 
wide class of DE models in order to constrain the dy- 
namics of Dark Energy form the observational data. As 
discussed in sec. I A we will propose an anzatz for L and 
y that covers the generic behavior of scalar field leading 
to an accelerating universe. 

For the quantity y = p m /V we propose to assume 
that V redshifts with an EOS w = w(z = 0) < 0, i.e. 
V = V (a /af { - 1+w °\ as in cq.(21). This does not mean 
that (x (a /a) 3 ( 1+w °> since the kinetic energy (p 2 /2, or 
cquivalently x, may grow faster or slower than V and the 
EOS for DE w(z) will be in general different than w - 
We have then 



V 



a 

Vo 1 — 



3w 



Vo{l + z) 



-3w 



(40) 



where p m = p mo (a /a) 3 , y Q = p mo /V = 20 m0 /f2^ o (l - 
w ) using V = P0o(l — w a )/2. Clearly the function y is 
an increasing function of z. From eq.(32) we know that 
as long as L is constant (or growing slower than y) that 
x and w will decrease as a function of z. 

If we want to have a constant EOS for DE at early 
times z > 1, as for example matter w = or radiation 
w = 1/3, which are reasonable behavior for particles, we 
should choose L proportional to y for large y or L/y — > 
if we want w — > — 1 at a large redshift. We then propose 
to take 



L = L a + Liytf(z) = L + L lVo 



1 + (a/a t ) k 



where we have chosen 



/(*) = 



1 



l + {a/a t ) k 1 + [(1 + *)/(! + *)]» 



(41) 



(42) 



and L Q , L±, z t , k are free constant parameters. The func- 
tion f(z) has as a limit f(z = 0) = [1 + (1 + z t ) fe ] _1 , 
f(z = z t ) = 1/2 and f(z ^> z t ) = 1. The parameters 
L a and L\ give w and w at large redshift z 3> z t while 
the transition epoch from w to w\ = w(z 3> z t ) is given 
by z t and k sets the steepness of the transition. The 
quantity £ takes the values £ = 1 and £ = only and we 
do not consider it as a free parameter but more as two 
different ansatze for L. 

Using eqs.(26) and (27) we can calculate easily x z = 
x y y z XLL z with 



L z = L y y z + Lff z = Li(y z f + yf z ), 
L y = £L 1 fi/~ 1 ,Lf = Liyt and for 



' 1 + *' 



fz = 



kp 



(1 + z) \ l + z 



Zt 



(43) 



(44) 



Eq.(43) becomes then 
(1 + *) 



L x = 



kf 



z t 



1 



- 3£w 



(45) 



with L\ f y£ = L — L a , or in terms of the scale factor a 
we have L a = (dz/da)L z = —L z (a /a) 2 and 



a/a \ \at 



3£,w 



(46) 



For £ = 1 eq.(41) allows a wide class of behaviors for w. If 
we want w to increase to w — 0, 1/3 we would take L\ = 
6, 12, respectively, or since in many scalar field models 
the evolution of w goes form w to a region dominated 
by the kinetic energy density with w — 1 and in this 
case we would should take L\ 3> 1. Of course a w(z ^> 
1) = 1 would only be valid for a limited period since fi^ 
should not dominated the universe at early times. We 
have included in eq.(41) the case £ = because we want 
to allow w to increase from w at small z and later go 
to w —> —1 (c.f. pink-dashed line in fig. (3)), since this is 
the behavior of potentials used as a models of DE as for 
example V = V 4>~ n ,n — 2/3 derived from gauge group 
dynamics [13] where the behavior of w(z) close to present 
time depends on the initial conditions. 



A. Initial Conditions and Free Parameters 

Let us summarize the parameters and initial conditions 
of our parametrization. The EOS w is only a function of 
x and a; is a function of L and y. From eq.(40) we see that 
y depends on two parameters w and 0,^, o (or equivalcntly 
on y , w ), since we are assuming a flat universe with DE 
and Matter and f2 TO0 = 1 — f2^ . From eqs.(41) and (40) 
we have the initial conditions at present time as 



Vo 



Pmo 



2 a 



(i 



L„ 



12(1+ wo) +6y {l-w 2 ) 

(I -Wo) 2 



l + (l + z t ) k 



(47) 



(48) 



For £ = 1, from eq.(35), the value of Lx gives the early 
time EOS, w t = w(z > z t ) = (L\ - 6)/(Zi + 6) or 
inverting this expression we have 



6(k;i + 1) 
1 — Wi 



(49) 



giving for example L\ = 6 for ioi = 0, L\ = 12 for 
W\ = 1/3, L\ =0 has w\ — > — 1 and Ii > 1 has w\ — > 
1. In the case £ = we have the limit L/y — >■ and 
w — > —1, independent on the values of L a , Li,z t , q. From 
eq.(41) we have that L depends on y and L Q , L\, z t , k and 
w , Qcf>o- However, not all parameters are independent, 
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FIG. 1: We show the evolution of fi^, to and for different 
models. We have taken £ = 1, k = 2,Li = 6 with 
zt = 0.1, 1,2, 10, 100 (red, dark blue, light blue, green and yellow, 
respectively). In black we have w,(l w using w in eq. (50) and Q C c 
for a cosmological constant (black dot-dashed). We take in all 
cases w = —0.9, f!^ D = 0.74 and L\ = 6 giving 
toi = w (z 2> zt) = for large z. 
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FIG . 2 : We show the evolution of fl* , w and for different 
models. We have taken £ = 1, zt = 2,Li = 6 fixed and k = 1/2, 2, 
5, 10, 20 (red, dark blue, light blue, green and yellow, respectively) 
In black we have w, Q w using to in cq.(50) and f! cc for a 
cosmological constant (black dot-dashed). We take in all cases 
w = —0.9, fl,/,,, = 0.74 and L\ = 6 giving w\ = w(z ^> zt) = for 

large z. 



since L a is a function of w , fl^o, L\,Zt,q and we are left 
with Vt^ and four parameters in w. To conclude, the 
free parameters are and for the EOS we can take w , 
wi, the transition redshift Zt and the steepness of the 
transition k. 



B. Other Parameterizations 

We present here some widely used parameterizations 
and we compare them with our DE present work. Let 
is present first a very simple and widely used DE 



parametrization [8] given in terms of only two param- 
eters 



w{a) + w a (l - a) = w Q + w a 



l + z 



(50) 



with a the derivative 



dw 
da 



dw dadw 

-w a , -j- = -T~T~ = v 1 + z ) w a (51) 
dz dz da 



Clearly w in cq.(50) is convenient since it is a simple 
EOS and it has only two parameters. However, it may 
be too restrictive and we do not sec a clear connection 
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FIG. 3: We show the evolution of Clj>, w and for different 
models. We have taken £ = 1, zt = 2 fixed and 
(k, Li) = (2, 6), (20, 12), (2, 100), (20, 1000), (20, 0) (red, dark blue, 
light blue, green and yellow, respectively) and £ = with 
z t = 0.1, k = 10, LI = 12 in pink-dashed line. In black we have 
using w in cq. (50) and Q cc for a cosmological constant 
(black dot-dashed). We take in all cases w = —0.9, f2<^ = 0.74 
and L\ = 6 giving w\ = w(z ^> zt) = for large z. 



between the value of w at small a and its derivative at 
present time w a . It has only 3 parameters fl^ 0l w ,w a , 
two less than our model but our model has a much richer 
structure. 

Another interesting parametrization was presented in 
[9]. It has 4 free parameters 



w = III 



+ (ioi - w a )G, G 



I _|_ e a d /d y — e {l-a)/di 



1 — e X l d 1 + e (a d ~a)l/d 

(52) 

where w a ,W2,ad,d are constant parameters. The func- 



tion G is constraint between < T < 1 with G = for 
a 3> ad and G = 1 for a <C a& Therefore ad is the scale 
factor where the transition of the EOS w goes from w to 
W\ The parameter d gives the width between the transi- 
tion, for small d the transition from w to wi is steeper. 
Even though w in eq.(52) gives a large variety of DE 
behavior [9], however the sign of the slope is fixed, and 
therefore our parametrization in eqs.(17) and (41) has a 
richer structure with the same number of parameters. 

In the interesting work of [10] they have followed a 
similar motivation as in the present work. They have 
presented a DE parametrization motivated by the dy- 
namics of a scalar field. Their parametrization has either 
two three parameters in eqs.(25) and (28) in the paper 
[10]. respectively (they do not take a eq as a free param- 
eter but we do think it is an extra parameter). The 
two parameter involves the quantities w;(a ^> a Q ), which 
gives the EOS at an early time, and X(a eq ) = V /V\ acq 
at DM-DE equality (i.e. f2 m = fide)- The second case, 
the parametrization also involves a term £ s (eq.(23) in 
[10]) which depends on seconde derivative of V and on 
the value of <j)/H at DM-DE equality. Since the func- 
tional form of the evolution of the EOS w(a/a eq ) is fixed 
in their parametrization the value of w at present time 
is determined if we know the value of a /a eq . Therefore, 
the quantity a eq must also be assumed as a free param- 
eter. As in our present work, they system of equations 
do not close without the knowledge of the complete V 
as a function of (j). However, since we are both inter- 
ested in extracting information from the observational 
data to determine the scalar potential the parametriza- 
tion given in eqs.(25) and (28) in [10] is a proposal to 
study a wide range of potentials V . Here we have taken 
a different parametrization which has a closer connection 
to the scalar potential V(4>) given by eqs.(9) and (8). 



C. Results 

We have plotted w for different sets of the parameters 
in figs. (1), (2) and (3) to show how w depends on L a , L\, 
Zt and k. We notice that our parametrization in cq.(41) 
has a very rich structure allowing for w to grow and or 
decrease at different rcdshifts. We have also plotted fi^ 
and the adiabatic sound speed c a defined in eq.(61) for 
each model. We are shoing some extreme cases which do 
not expect to be observational valid but we want to show 
the full extend of our parametrization. We have also in- 
clude a cosmological constant "CO" (black dot-dashed) 
in the figures of fl^ and since w cc = — 1 and c acc = 
we do not include them in the graphs for w, (? a . We also 
plotted Wjfl^^c^ (in black) for the parametrization in 
eq.(50) for comparison. We take in all cases w Q = —0.9, 
= 0.74. 

In fig.(l) we show the evolution of 0,$, w and c a for 
different models. We have taken £ = l,k = 2, L\ = 
6 with z t = 0.1,1,2,10,100 (red, dark blue, light blue, 
green and yellow, respectively). We take in all cases L\ — 
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6 giving wi = w(z z t ) — 0. Notice that the yellow line 
the increase to w = does not show in the graph since 
Zt = 100 and we plot w only up to z = 4. The slope in 
w depends on the value of zt and since k = 2 is not large 
the transition is not too steep. The value of fi^ decreases 
slower than a C.C. and for smaller zt it decreases even 
more slowly, i.e when w approaches zero faster. 

In fig. (2) we show the evolution of f2^, w and c 2 
for different models. In this case we have taken £ = 
l,z t = 2, L\ = 6 fixed with w± = and we vary 
k = 1/2, 2, 5, 10, 20 (red, dark blue, light blue, green and 
yellow, respectively). We clearly see in the evolution w 
how the steepness of the transition depends k and that 
c 2 has a bump at zt and it is more prominent for steeper 
transition. This is generic behavior and we could expect 
to see a signature of the transition in large scale struc- 
ture. 

In fig. (3) we show the evolution of fi^, w and c 2 
and we take £ = l,Zt = 2 fixed and (k, L\) = 
(2, 6), (20, 12), (2, 100), (20, 1000), (20, 0) (red, dark blue, 
light blue, green and yellow, respectively) and £ = with 
Zt = 0.1, q = 10, LI = 12 (pink-dashed line). In this case 
we vary L\ and we see that for large L\ the EOS w be- 
comes bigger and it may approach w ~ 1 (e.g. green 
line). Of course this case is not phenomenologically vi- 
able but we plot it to show the distinctive cases of our 
w parametrization. Once again, a steep transition gives 
a bump in c 2 . The pink-dashed line shows how w can 
increase at low z and than approach w = — 1. 

We have seen that a our parametrization gives a wide 
class of w behavior, with increasing and decreasing w. 
From the observational date we should be able to fix the 
parameters of L in eq.(41) and we could then have a 
much better understanding on the underlying potential 
V{4>) using L = (V'/V) 2 . 



IV. PERTURBATIONS 



Eq.(53) can be expressed as a function of a with Y = 
aHY a for Y = 5<p,<p,H,h and the subscript a means 
derivative w.r.t. a (i.e. Y a = dY/da), giving 



5<t> = --jh a (j) a . 









' k 2 V"' 


(1 


" H J 


_a 2 H 2 + H 2 



In the slow roll approximation we have 
V" 

Tx < 3 



3H 2 

where we have used eq.(17) and 

V"V 



r = 



yr. 



(54) 
(55) 

(56) 



Eq.(55) implies that an EOS of DE between — 1 < w < 
-1/3,0, 1/3, with < x < 1/2,1,2, requires r < 3/x — 
6,3,3/2, respectively. For a scalar field <f> to be in the 
tracking regime one requires Y to be approximated con- 
stant with r > 1 [11]. Therefore the regime 1 < Y < 3/x 
allows a tracking behavior satisfying also the slow roll ap- 
proximation. Here we are more interested in the late time 
evolution of DE and the tracking regime is not required 
and in fact we expect deviations from it. However, if Y is 
nearly constant the evolution of the perturbations in (54) 
are then given only in terms of x and we can use our DE 
parametrization in eq.(41) and eq.(56) and to calculate 
them. 

We can express the slow roll parameter e, T in terms 
of r and L as 



L 
2' 



V" 

t = - = yl. 



(57) 



We have decided to use L, T instead of e, r) not to confuse 
the reader with the inflation parameters and the DE ones. 



Besides the evolution of the homogenous part of Dark 
Energy </>(£), its perturbations 5(p(t,x) are also an es- 
sential ingredient in determining the nature of DE. The 
formalism we work is the synchronous gauge and the lin- 
ear perturbations have a line element ds 2 = a 2 (—dT 2 + 
(Sij + hij)dx t dx^ , where h is the trace of the metric 
perturbations [17, 19]. In this sect. (IV) a dot repre- 
sents derivative with respect to conformal time r and 
H = a/a = (da/dr)/a is the Hubble constant w.r.t. t, 
while H = (da/dt)/a 



Fluid Perturbations 



The evolution of the energy density perturbation 5 
Sp/p, 8 the velocity perturbation [17-20] 



-3H(cj-iu 



S = -(l + w )(k 2 + 9H 2 [c 2 ~c 2 a })^ 

S 



1 + w 



h 
2 



(58) 



Scalar Field Perturbations 



2 ;„2 



9 = -H(l - 3cJ)0 + cik 



1 + w' 



(59) 



For a DE given in terms of a scalar field, the evolution 
requires the knowledge of V and V while the evolution 
of S(j)(t, x) needs V" through [17, 19] 



+ 2H(50 + [fc 2 + a 2 V"\ 



1 ,■ 



-h(j). 



(53) 



and we do not consider an anisotropic stress. The evo- 
lution of the perturbations depend on three quantities 
[19, 20] 



P 



(60) 
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P 
P 



5p 



w + w— 

p 

w 

3H(1 + w) 



= w + 



x z w x 
3a(l + w) 



(61) 



(62) 



where w is the EOS, H the Hubble constant in conformal 
time, c 2 is the adiabatic sound speed and c 2 is the sound 
speed in the rest frame of the fluid [18, 19]. For a perfect 



fluid one has 



c 2 but scalar fields are not perfect 



fluids. The entropy perturbation Gi for a fluid Pi with 

Si = Sp./pi are 

WiGi = (c si - c ai )Si = — (- -) (63) 

Pi Pi Pi 

where the quantities Gi and cr ai are scale independent 
and gauge invariant but (? si can be neither [17, 20]. In its 
rest frame a scalar filed <j) with a canonical kinetic term 
one has c 2 = Sp/Sp = 1 [18, 19]. One can relate the rest 
frame 6,9 to an arbitrary frame S, 8 by [20] 



6 = S + m(l + w) 



k' 2 



(64) 



and 



Sp = c 2 8p + {c 2 



cl)m{l + w) Pl - (65) 



As we see from eqs.(58)-(59) the evolution of 8 depends 
on c 2 , c 2 and w. Using eq.(60) and since w is a function 
of x(a) we have 

2 . X Z W X 

c a = w + TTT, — T 66 
3a(l + w) 

with w x /(l + w) = 2/[x(l + x)]. From eqs.(44) and (45) 
we can express c 2 as a function of the parameters of x. 

Finally, we can relate V" in terms of the adiabatic 
sound speed c 2 in cq.(61) and its time derivative, using 
c\ = p/p = 1 + 2V'/3H((fy/dt), giving 



dt 




\3H 2 



and for c 2 ^ 1 we can invert eq.(67) to give 
V" TV _ 1 dc 2 a , 



(67) 



(68) 



where we have used H = dH/dt = —(pr +Pt)/% = 
—3H 2 (1+wt)/2 with pt,Pt, wt the total energy density, 
pressure and EOS, respectively. In our case we have P t = 
Pm + P4>, PT — Pm + P<j> — P(j> and using p^ = V(l + x) 
and eqs.(8) and (12) we have 



_ Pt n wp^ 

Wt — = Wild, = 



Pt 



3H 2 



w(l + x) 
1 + x + y 



1 



1 



+ y 

(69) 

With eq.(69) the l.h.s. of eq.(68) depends then only on 
y, x and are fully determined by our parametrization. In 
the full slow roll approximation = and one has c 2 a = 
-1. 



V. CONCLUSIONS 

We have presented a new parametrization of Dark En- 
ergy This parametrization has a rich structure and al- 
lows for w(z) to have a wide class of behavior, it may 
grow and later decrease or other way around. The 
parametrization of w is given in terms of x(L, y), given in 
eqs.(27), (41) and (40). The EOS w is constraint between 
— l<w<\ for any value of x, with < x by definition. 
The free parameters of w are L a , L\, z tl k, or alternatively 
w and the EOS at an early time w\ = w(z 3> z t ), given 
by L a and L%, respectively (c.f. cq.(48), while z t gives 
the transition redshift between w and w\ and k sets the 
steepness of the transition. Besides studying the evolu- 
tion of Dark Energy we also determined its perturbations 
from the adiabatic sound speed c 2 a and c 2 given in eqs.(61) 
and (62), which arc functions of x and its derivatives. We 
have seen that a steep transition has a bump in c 2 a and 
this should be detectable in large scale structure. 

We can use the parametrization of x(L, y) in eqs.(27), 
(41) and (40) and c 2 , and c 2 in cqs.(61) and (62) without 
any reference to the underlying physics, namely the dy- 
namics of the scalar field (j>, and the parametrization is 
well defined. However, it is when we interpret x = 4> 2 /2V 
and L = (V'/V) 2 A and y — p m /V that we are analyz- 
ing the evolution of a scalar field cf> and we can connect 
the evolution of w to the potential V((f>), once the free 
parameters are phcnomcnological determined by the cos- 
mological data (only when we take \q\ 3> 1, A ~ 1 are we 
taking in the slow roll approximation). 

To conclude, we have proposed a new parametrization 
of DE which has a rich structure, and the determination 
of its parameters will help us to understand the nature 
of Dark Energy 



ACKNOWLEDGMENTS 

We acknowledge financial support from Conacyt 
Proyect 80519. 

Appendix A 

The parameter \q = 4>/V'\ is clearly smaller than one 
in the slow roll regime (</> < 3H<fi ~ V'), and let us now 
determine the dependence of q on the potential V(<p) and 
its derivatives. The evolution of q is 



1 yi y,2 



3H 



H 



-q+(l + q)T7T^ + 2Tx ) 



where we used eq.(14), 

'4> u'> 3H<j) m'cj) 



V 



V 



V 



V 



(Al) 
(A2) 

(A3) 
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-3Hq + 3H(l + q) 



H 



9H 2 3H 2 



(A4) 



and 



V" 
9H 2 



Tx 



V"V 

yi2 ■ 



(A5) 



We can estimate the value of q if we drop the term pro- 
portional to <fi in cq.(A3) and using eq.(37) giving 



V" 



H 
IH 2 



1 - f-^T 



H 
IIP 



(A6) 



In a stable evolution of (j> we have a positive V" and 
since H is negative both terms have opposite signs, but 
of course we do not expect a complete cancelation of these 
terms. However both of them are smaller than one, since 
< -H/3H 2 < 1/2 for x < 1 and \V" /9H 2 \ = Tx/3 < 
1/3 in the slow roll approximation. A tracker behavior 
requires T > 1 [11] and x < 1/T < 1. Finally, the evo- 
lution of L is given by L = 2AA(1 + q) 2 + A 2 g(l + q) = 



2L[\/\ + q(l + q)}, 
L = — y — + 



2Lq 



= 3HL 



2x 



3HL 



(1 + q) (l + q) \l+q 

2{q - x){l + 2x + 2y) - y(l + q) 



H 
3H~ 2 



2(1 + <?)(! + x + y) 



(A7) 



We see that at — 1 < q < x we have L < giving a 
decreasing L as a function of time or an increasing L as 
a function of z. For (q — x) / (1 + q) > y/(l + 2x + 2y) or 
equivalently for q > (y + 2.x (1 + 2x + 2y))/(2 + 4x + 3y) 
we have L > and a decreasing L as a function of z. The 
evolution of d(L/y)/dt is given in eq.(36). 

Instead of choosing a DE parametrization as in eq.(41) 
we could solve eqs.(Al) and (A7) for different potentials 
V{4>) or by taking different approximated solutions or 
ansatze for T. However, we choose to parameterize di- 
rectly L as in eq.(41). Still using L = aL a = aHL a and 
from eqs. (46), (A7) and L^f = L - L Q we identify 



(L-L )\kf [-) -3£w 



3L 



2x~q H 



l + q 3H 2 I 

(A8) 

and the choices of T and q would fix the parameters L Q , k 
and at- 
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